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An integrable deformation of the known integrable model of two interacting p-dimensional and 
g-dimensional spherical tops is considered. After reduction this system gives rise to the generalized 
Lagrange and the Kowalevski tops. The corresponding Lax matrices and classical r-matrices are 
calculated. 
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^ ! 1 Introduction 

. ■ 

The most common examples of classical .R-matrices are associated with decompositions of 
Lie algebra into a direct sum of two Lie subalgebras £Q . According to 2\ we can consider 
various "perturbations" of the standard decompositions of the loop algebras which may be 
associated with integrable deformations of the known integrable systems (HI E] • In this note 
we consider integrable deformations of some tops associated with the Lie algebra so(p, q) 
and calculate the corresponding /^-matrices. 

Let be a self-dual Lie algebra, 0+,0- C g its two Lie subalgebras such that g = 
+ -i-g_ as a linear space. Let P + , P be the projection operators onto g± parallel to the 
^ ! complimentary subalgebra; the operator 

/? P + -P_ (1.1) 

^ ■ is a classical i?-matrix. If a is an intertwining operator 

^ ; a[X, Y] = [aX, Y] = [X, aY] (1.2) 

for any X, Y G Q, then R a = R o a is also a classical i?-matrix jT]. So, one decomposition of 
determines a family of /^-brackets those orbits are related 

L^a^L, (1.3) 

if a is invertible. For instance, such change of variables was used in jT] to connect Lax 
matrices for the Manakov top on so(4) and for the Clebsch system on e(3). 
The most interesting class of examples is provided by loop algebras 5£ 

££ = q[\, A" 1 ] = { X(X) = XiEa\. (1.4) 
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The standard decomposition Jz? = Jzf + -j-Jzf_ is defined by the natural Z-grading in powers 
of auxiliare variable A (spectral parameter) 

^+ = 00*^, if- = 0ftA' (1.5) 

i>0 i<0 

The i?-bracket associated with decomposition (|1.5|) has a large collection of finite-dimensional 
Poisson subspaces 

They are ae^-invariant subspaces 

ad* R X ■ L = (ad* s X + ■ L)+ - (ad*X_ • L)_, X± = P±X, L G J£ mn , (1.6) 

which are invariant with respect to the Lax equation 

L t = [L,M], M = P + (dH(L)). (1.7) 

Intertwining operators in Jzf are multiplication operators by scalar Laurent polynomials. 
The crucial observation in [3J 0] is that to construct new Lax matrices by the rule ()1.3j) we 
can use the matrix Laurent polynomials a. 

According to [2], if Q is associative algebra, the same space 5£ (jl.4j) is a linear sum of 
(|1.5|) and subalgebra JzfZ 



I i=i 



li 



defined by some fixed element y G g. This new decomposition of JSf defines another classical 
i?-matrix (fTTj) 

Ry = p y_py. (l.g) 

Subalgebra j£ff could be regarded as some "perturbation" of the standard subalgebra J??_ 
by element y. In this approach the finite-dimensional Poisson subspaces of the ^-bracket 
(JIl) 

n n 

, A* (1 - X^y) or (1 - X~ 1 y)~ 1 g, A\ 

i=—m i=—m 

are deformations of the known orbits S£ mn of i?-bracket, which describe integrable deforma- 
tions of the known integrable systems. 

Our aim here is to show mapping of the Lax matrices 

L = Xa + l + X^s, M = Xb (1.10) 

associated with Cartan type decomposition of g = so(p, q) into the new Lax matrices asso- 
ciated with i? y -bracket 

L v = (l-\- 1 y)- 1 (\a + l' + \- 1 a), M y = M(l — X~ x y), (1.11) 

which describe new integrable deformations of so(p, q) tops. 
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2 Interacting spherical tops on so(p, q) algebra 

Let G = SO(p,q) is the group of pseudo-orthogonal matrices with signature (p, q), p > q. 
The Lie algebra so(p, q) consists of all the (p + q) x {p + q) matrices satisfying 

X T = -JXJ, 

where J = diag (1, . . . , 1; —1, . . . , — 1), tr J = p — q and T means a matrix transposition. 
In the natural (p, g)-block notation an element X G so(p, q) has the form 

X=( i S ), (2.1) 
\ s m J 

where l T = —£, m T = —m are p x p and q x q matrices, s is an arbitrary p x q matrix . 

The Cartan involution is given by aX = —X T and g = f + p is the corresponding Cartan 
decomposition. The maximal compact subalgebra f = so(p) © so(q) consists of matrices X 
with s = 

f = S0(p) © S0(q) = { ( I ° )} 
The subspace p consists of matrices X with I = 0, m = 




The pairing between so(p, q) and so(p, q)* is given by invariant inner product (X, Y) = 
— |trXY which is positive definite on f. 

We extend the involution a to the loop algebra Jf(g, a) by setting (aX)(X) = a(X(—X)). 
By definition, the twisted loop algebra Jif(so(p, q). a) consists of matrices X(X) such that 

X(\) = -X T (-X). (2.2) 

The pairing between Jz? (g, a) and (q, a) is given by < X, Y >= Res A _1 (X, Y). 



2.1 Quadratic Hamiltonian 

Let k G SO(p), r G SO(q) and £ 6 so(p), m G so(g) be the configuration and momentum 
variables on the phase space T*SO(p) x T*SO(q). All the non-zero Lie-Poisson brackets are 
equal to 

.. P \ — A. /. + A. />. _ A. p. _A. /. 
• • t 1 A- h ■ — /> t • 

(2.3) 

, fn^ n } Ttljfc -\- 5jk 'Trim ^jfc TTbjn ^j'n ^ifc 

The Poisson mapping from T*SO(p) xT*SO(q) into ^f(so(p, g), a) is defined by the following 
theorem pQ. 
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Theorem 1 (Reyman, Semenov-Tian-Shansky) 
The spectral invariants of the Lax matrix 



Xa + 1 + X x s 



kAr T 
rA T k T 



(2.4) 



A - 



klk 1 




rmr 



+ A- 1 



F 
F T 



are in involution with respect to the canonical Poisson brackets \2.'J\) on T*SO(p) x T*SO(q). 
Here A and F are constant p x q matrices. 

We refer the reader to \X\ for a complete proof of this theorem which uses that L(X) (J2.4|) is 
an orbit of classical P-matrix associated with a standard decomposition of the twisted 
loop algebra Jzf (so(p, q),cr). 

Here we present an elementary proof using a tensor form of the P-bracket which is more 
familiar in the inverse scattering method jT]. According to jH] functions <f>k{L) = tr L h , k > 2 
are in the involution if and only if matrix L(X) satisfies relation 



L(A),L(/i) 



r 12 {X,fj),L(X) 



r 2 i(A,/x) ,L{n) 



(2.5) 



Here we used the standard notations 



L(A) = L(A)®1, L(/j,) = 1®L(ij,), 

and matrices r 12 , r 2 i are kernels of the operators R and R* such that 

r 2 i(X,fi) = nri 2 (/i, A) n, 

where n is a permutation operator UX (g> Y = Y <8> XU. 

One checks without difficulty that L(X) ()2.4|) satisfies (j2.5J) with the following kernel of 
the i?-matrix 

r 12 (A, ^ = -^jPf + p^n, (2-6) 

where and P p are kernels of the standard Casimir operators acting in the orthogonal 
subspaces f and p respectively. They are equal to 



P (p-l) + q(q-l) 
2 



P'J 



a=l /3=1 

where antisymmetric Z a and symmetric Sp matrices form two orthonormal basises in f and 
p. For instance we can use 
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Corollary 1 The equations of motion generated by the spectral invariants of the matrix L(\) 
\2.4\j are Hamiltonian equations with respect to canonical brackets which give rise to 

Lax equation 7| ). 

Let A be the truncated diagonal matrix A = aidij. The spectral invariants <pk{L) of L(X) 
f!2.4|) may be combined in the following Hamiltonian 



a.jbj — a~bi 



cq - a"- 



V 



af -aq 



(2.7) 



i=g+l j=l 



which formally describes an interaction of p-dimensional and q- dimensional spherical tops. 
Here a^, 6j, i = 1 . . . q and 7 are arbitrary parameters jT]. In this case the second Lax matrix 
M in (jl.7j) is equal to 



M(A) = P+(dH(L)) = - 
where diagonal matrix B = biSij includes parameters \. 



kBr T 
rB T k T 



A, 



(2- 



2.2 Integrable deformations 

Let us introduce matrix 

y = c ( _° pT cec, (2.9) 

which is a generic solution of the following equations 

yXs = a(yXs), VlGfl, s = ( A ) G P ' (2 ' 10) 

This solution is parametrized by an arbitrary numerical parameter c. 

The matrix y (j2.9j) . y G so(p + 5) = f + ip, does not belong to the algebra so(j>, q) (|2.1)l . 
So, to consider deformations j£f^ (jl.8j) of J??_ we have to embed initial algebra q and the 
element y into some algebra and only then to discuss mapping of the known orbit (jl.lO|) 
of ^-bracket (H3J) into the orbit flTTTl) of i? y -bracket (Q defined in Jzf(fl)*. 

Therefore, for brevity, we shall consider kernels r(A ;/ u) and r y (X,fj,) instead of the cor- 
responding operators R (jl.lll.5j) and R y (jl.9ll.8j) . 
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Proposition 1 The spectral invariants of the Lax matrix 

L y = (1 - \- l y)-\\a + V + A _1 s), 

where 



l ' = 1 ~ \i ya + ay \ 



(2-11) 



f2.12l 



are m involution with respect to the canonical Poisson brackets \2.'J\) on T*SO(p) x T*SO(q). 
To proof this proposition we have to check equation (J2.5j) with the kernel of the operator R y 



rf 2 (A, //) = [l ® (1 - /T^)" 1 ] r 12 (A, /z) [(1 - A y) ® 1 
by using properties of the element y (I2.1(J|) and definition of matrices ty, (A, /i) (|2.6j) . 



(2-13) 



Corollary 2 T/ie equations of motion generated by the spectral invariants of the matrix L y 
H2.ll)) are Hamiltonian equations with respect to canonical brackets h2.ty . They give rise to 
Lax equation where 



M y = Pl{dH{L)) = M{l-\- 1 y). 
The Lax equation (|1.7|) for L y (j2.11j) may be rewritten as a Lax triad 

^-L l = L i ML 2 -L 2 ML u 2 = 1,2 
at 

on the pair of matrices 

and L 2 = 1 — A -1 ?/, 



(2-14) 



(2-15) 



Lx = L- -(ya + ay) 

entering in the definition of L y ([2. lip . 

For future reference we present two equivalent forms H± 2 of the deformed Hamiltonian 
H y using additional canonical transformations of the phase space 



I — > I ± - (sa — as). 

2 v ; 



(2.16) 



Compositions of these canonical transformations with noncanonical map (|2.12j) 

I _> l 12 = I _ I zp cs )a + a(y ± ca)) 
act either in the subalgebra so(p) 

£ _, ^ = £ - c (k T F r A T - Ar T F T k s j , 
either in the subalgebra so(q) 

m — > mi = m — cyr A k F — F k Ar 



(2-17) 



f2.18) 



The corresponding perturbations of the initial Hamiltonian H (J2.7|) depend either on 
variables 



H x = H -c-ti(k£Br T F T ) + —ti(B T ArF T Fr T ) 



either on m variables 



(2.19) 



H 2 = H + c ■ tr (r m B 1 k 1 F) + — tr (B A 1 k F F k) . 



The corresponding Lax matrices are equal to 



or 



4 1} = a - A" 1 *)" 1 



4 2) = (i - ^yr 1 



L(X) + c 



L(X) - c 



FrA T k T - kAr T F T 



















FrA T k T -kAr T F T 



(2.20) 



(2.21) 



(2.22) 



We have discussed so far Lax matrices for the tops in the stationary reference frame, 
i.e. used Euler-Lagrange description of motion. According to P, we can go over to the Lax 
matrices in the frame moving with the body which amounts to the gauge transformation 



L(X) 



g'Lg = \h + / + A s 



Here g 



A 
A T 



£ \ j / k T Ft 
m + [ r T F T k 







(2.23) 



M(A) = g T M(A)g + w = -[ B ° T q]a + g T j t g. 



k 
r 

m by a linear transformations. 



and w = g T g is the block-matrix of angular velocities, related to £ and 



Proposition 2 The Lax matrix L(X) associated with the Euler-Poison description of 

the motion satisfies equation \2. 5}) with the following matrices rij(X,fi) 



ri 2 (A,/x) = -r 12 (A : ), r 2 i(A,/x) = -r 2 i(A i ), 
which are obtained from r 12 (A,/x) \2. 6}) by change of the spectral parameters. 



(2.24) 



Using the similar gauge transformation to the matrices L y (X) (J2.ll)) and M y (X) (|2.14j) 
one gets 



1 



L y (X) = g 1 L y g={l-X- 1 y)- 1 [L{X)--{y~a+~ay) 



M y (X) = g T M y (X)g + g T j t g = M(l-X- 1 y). 



(2.25) 
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Here matrices w y = g T g are calculated with respect to cubic Hamiltonian H y and, therefore, 
it depends on all the dynamical variables. As above (|2.15|) . a pair of matrices 

L\ = L + -{ya + ay) and L 2 = 1 — A _1 y, 

entering in the definition L y ()2.25j) . satisfy to the Lax triad 

^Uf(A) = Li(A) M(A) + M T (— A) L,(A), z = 1, 2. (2.26) 
at 

These equations ()2.26|) have a more complicated structure with respect to Lax triad in the 
rest frame (j2.15|) 

In the body frame the matrix 



k T Fr 



» g S I _ r T F T k 

depends on the dynamical variables and, therefore, the kernel (j2.13j) of the corresponding 
operator R y depends on these variables too. 

However in the moving frame we can construct another Lax matrix such that the kernel of 
the corresponding i?-matrix does not depend on dynamical variables. Namely, using matrix 

5 = c ( _° A T A ) (2.27) 

such that 

zXa = a(zXa), VAT E g, a = ( ® T ^ J G p, (2.28) 
we can consider perturbation of the Lax matrix L(X) (|2.2Hj) by z 

L Z (X) = (1 - A5)- 1 (Z{\) -hzs + sz)). (2.29) 



Proposition 3 The Lax matrix L Z (X) satisfies equation h2.ty with the following nu- 

merical r-matrix 



r 1 z 2 (A, / u)= l^l-AFr 1 r 12 (A,//) (1 - \iz) <g> 1 



(2.30) 



The proof is straightforward. 

Applying inverse gauge transformation to L 2 (A) one gets another Lax matrix L Z (X) in 
the stationary frame 



L z (X)=gL z g 7 



'1 - Xz) 



1 (l(\)-±(Z8 + 8z)\ 



z = gzg T . 



(2.3L 



The corresponding r-matrix is dynamical. 

So, for initial Hamiltonian H (|2.7j) we know two Lax matrices L(X) ()2.4j) and L(X) 
()2.23|) associated with the physically different coordinate systems. For these Lax matrix we 
constructed by two perturbations L y (j2.11|) . L z ()2.31|) in the rest frame and L y ()2.25|) . L z 
(12.29)1 in the body frame. 



3 Further examples 



According to pQ , we shall try to exclude the nonphysical degrees of freedom using symmetry 
of the Hamiltonian H (12. 7j) and its perturbations Hi j2 (j2.19l2.2(Jj) . Assume that A = E is 
the truncated identity matrix E^ = 5ij. 

Below all the Hamiltonians will be expressed through kinetic momentum £y G so(p) and 
the entries of the Poisson vectors Xi = k T fi, where fi are the column vectors of the matrix 
F. These variables are canonical coordinates on the algebra e(p, q)* = so(p) x {BP ® W 1 ) pp. 

3.1 Algebra so(p, 1), Lagrange top and spherical pendulum. 

If q — 1 the phase space is T*SO{p). We can put 

m = 0, r = l, (3.1) 

without loss of generality. In this particular case the deformations (j2.19j) and ()2.2()j) of the 
initial Hamiltonian are quadratic polynomials instead of cubic ones. 
Following to 0, we can rewrite Hamiltonian ()2.7)1 in the form 

tf = ^E4 + iE4-( e ^)> x = kT f- ( 3 - 2 ) 

i,j=l i,j=2 

Here e\ is a first vector of the standard basis in W. 

The Hamiltonian ()3.2j) describes rotation of a rigid body around a fixed point in a 
homogeneous gravity field. The vector x is the vector along the gravity field, with respect 
to the body frame and e\ is the vector pointing from the fixed point to the center of mass 
of the body. 

It is Lagrange case because the body is rotationally symmetric and the fixed point lies 
on the symmetry axis e\. The perturbations (j2.19|) and ()2.20|) of the Hamiltonian (j3.2j) 

2 2 

H l = H+ C -{e u k T ff = H+°- (e l5 xf 

H 2 = H + c(£ei, k T f) = H + c(£e h x). 

are quadratic polynomials and, therefore, the corresponding equations of motion have the 
form of the Kirchhoff equations. 

Let x = k T f and 7r are canonical coordinates in M? p , {iTi,Xj} = 5ij. Substituting = 
XiiTj — TCiXj in L(\) ()2.4j) one gets a Lax matrix for the spherical pendulum in appropriate 
physical coordinates pQ. The corresponding Hamiltonian 

describes a motion on the sphere (x, x) — 1, (x, ir) = 0, in a homogeneous gravity field. 

Its deformations flZZD} and fl2ZG2D look like 

H 1 = H + c^Aiin, H 2 = H+ °— A i x i) ■ 
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3.2 The Kowalevski top. 

According [TJ |3] , the generalized p-dimensional Kowalevski top is the reduced system (|2.7|) 
with respect to the action of the subgroup SO(q). The reduction amounts to imposing the 
constraints 

m + P£P = 0, r = l. (3.3) 

Here P means the orthogonal projection from W onto subspace of R 9 spanned by first q 
vectors of the standard basis. The reduced phase space is T*SO(p) with its canonical Poisson 
structure. 

Inserting the constraints (j3.3|) into ()2.4|) one gets the Lax matrix for the generalized 
Kowalevski top 

tKow/\\ _ ( kE\ f k£k T \ +x -i( F\ 

L {X) -{E T k T ) X ~{ -E T tE) +X \F T Oj- (3 - 4) 

Proposition 4 The spectral invariants of the Lax matrices L Kow (X) \3.4\) and L y (X) \2.11\) 
by \3. 3]) are in the involution with respect to Lie-Poisson brackets on T*SO(p). 

We present a direct proof of the involutivity of invariants without recourse to Hamiltonian 
reduction. 

In the rest frame it can easily be shown that the reduced Lax matrices satisfy equation 
()2.5|) with the same kernels ri 2 (A,/i) ()2.6|) and rf 2 (A,/i) f)2.13j) . These kernels do not change 
by the Poisson reduction (|3.3|) and the corresponding operators R (jl.lj) and R y (|1.9|) remain 
differences of the same projectors. 

It should be pointed out that constraints ()3.3|) agree with the Euler-Lagrange description 
of the top, but disagree with its Euler- Poisson description pQ. In the body frame the Lax 
matrix is given by 



Proposition 5 The spectral invariants of the Lax matrices L (X) \3. 5\) and L Z (X) \2.i 
by \3. 3]) are in the involution with respect to Lie-Poisson brackets on T*SO(p). 

In the body frame reduction ()3.3|) does not Poisson mapping which changes /^-brackets. 
Thus, after reduction the Lax matrix L(X) (j2.23|) satisfies (|2.5J) with the reduced kernel 



p(p-l) + qr(qr-l) 
2 



f^(A,/i)=r 12 (A,/i)- J2 Z a ®PZ a P, 



(3.6) 



a=l 



Here 



P 



E 
E T 



PZ a P = { 



Z * + 
0. 



9(1,-1) 



1) , < a < 



2 ' 



p(p-i) 



< a. 
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Algebro-geometric description of the corresponding "non-standard" R operator may found 
in0. 

After reduction ()3.3j) the perturbed Lax matrix L 2 (A) ()2.29j) satisfies ()2.5|) with the 
reduced r-matrix 

(1 - fiz) ® 1 . 

Substituting constraints (j3.3|) into ()2.7|) gives the reduced Hamiltonian on T*SO{p) 

\i,j=l «,i=l / i=l 

where ex, • . . , e p is the moving orthonormal frame associated with the top, and fx, . . . , f q are 
the column vectors of F. This Hamiltonian describes a p-dimensional top with a very special 
SO(q) x SO (p— ^-symmetric inertial tensor, rotating under the influence of g constant forces 

P- 

Inserting constraints ()3.3|) into ifi ()2.19|) and #2 (|2.2U|) one gets quadratic polynomials 

Hl = H-c J2(P£P ei , k T f,) + ^ X> T /,, fc T /,) , (3.8) 
i=i i=i 

and 

H 2 = H-c^^k T f l ), (3.9) 
i=i 

It means that the corresponding equations of motion have the form of the Kirchhoff equations 
of motion of a rigid body in the ideal fluid. 

Thus, the characteristic doubling of the terms £ 2 -, for i,j^q in (13 .7j) . and existence of 
the nontrivial quadratic perturbations Hf (|3.8|) and H 2 V ()3.9|) are results of reduction with 
respect to SO(q). 

The case p = 3, q = 2 and /2 = gives usual Kowalevski top on the algebra e(3) with a 
Poisson vector x = k T fx. The deformations (|2. 19112.21111 of the standard Hamiltonian became 

1 / s c 2 

#i = 2(4 + 4 + 2 4j -x l + c£ l2 x 2 + -{x\ + x 2 2 ). (3.10) 

and 

#2 = ^(4 + 4 + 24) -xi + cfe^ + teO- (3.11) 

Recall, ifi and H 2 are two different forms of the one Hamiltonian H y , which are connected 
by canonical transformation ()2.16|) . 

3.3 Integrable systems on so(4). 

If p = q = 3 and Fy = the Hamiltonian H (|2.7|) describes integrable system on the phase 
space so(4) = so(3) © so(3). Let 

^fe ^ijk^-jki ^k ^■ijk'^jki 



zKow 



12 



(A,//) 



+ - A5i 



12 



(A,//) 
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where is the totally skew-symmetric tensor. The Hamiltonian (|2.7|) becomes 



UkVk- 



The coefficients 



Oik 



k=l 



2 CLj^b>i ^j^j 



k=l 



(3 k 



2 (Xjjjj Qjj 



satisfy the relation 

(«i - a 2 )0l + (a 2 - a 3 )/3f + (a 3 - a^ffi + {a x - a 2 )(a 2 - a 3 )(a 3 - a x ) = 0. 

Thus, according to jS], this integrable system belongs to the so-called Steklov-Manakov 
family of integrable systems, characterized by the property that there exists an additional 
quadratic integral. In this case F = and proposed deformations (|2.19H2.2T7|) are trivial. 
The more complicated passage to the algebra so(4) was proposed in jlj. Let us apply 

(o\ ______ 

composition of the gauge transformations to the Lax matrix L y ()2.22j) 

L y = (l-\- 1 9 )(g T LV(\) g )(l-\-i g )-\ g = c(° Q 

Substituting constraints (|3.3|) one gets 

-i -l 



1 + 



A 2 




F T F 



x 



(3.12) 



E 
E T 



A - 



f 







—E T £ E 



+ A- 1 



where 
The matrix 



h 




F T F 



^F 
F T 1 



1= (l + f)k and d = k T £k. 





r T F T kk T Fr 



c 2 




F T dF 



(3.13) 



is numerical because r = 1 ()3.3j) only. 

Proposition 6 The Lax matrix L y {\) satisfies \2. 5)) with numerical r -matrix 



1® ( 1 + J l 2 h 



-1 



l +Y2 h 



(3.14) 



The proof is straightforward. Thus, for the generalized Kowalevski top we constructed third 
Lax matrix with the numerical r-matrix. In contrast with the previous two Lax matrices in 
this case we do not know the physical meaning of the corresponding coordinate system. 
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The non-zero Poisson brackets between variables £, 7 and d (|3.13j) are 

.. p \ = 5. p. 4. s £ -5- £■ - 8- £■ 

ij ; ^nm } ^j'm 7m ^im Tnj ) 

(ifc n } = dj n + 5j n (ijfc 5j n djfc d)^/. dj n (3.15) 

\ d%j ; 7«fe J ^ira Ikj &jn ^ki > 

{7«j;7fcn} ^ {.&jn,£is &ikdjn) • 

Entries of matrix £ jointly with entries of any column of 7 form subalgebras with respect to 
bracket (J3.15|) . If rankF = 1, then we can put 

F T dF = 0, ( 7 T F)y = 5 jllu 

and, therefore, the Lax matrix ()3.12|) is defined on the one of such subalgebras only. 
Let p = 3, in canonical variables 

Ji &ijk£jk ; Uj 7jl 

the brackets ()3. 1 5|) coincide with the standard Lie-Poisson brackets on so(4) 

{Ji,Jj} = EijkJk , { Ji , yj } = e ijk y k , {yi ,yj} = c 2 e ijk J k . (3.16) 

Thus, subalgebra {£, 7,1} is isomorphic to the Lie algebra so(4) [I]. 
If q = 1 or g = 3 the perturbed Hamiltonian ()2.20|) 

H 2 = Jl + 4 + J*- yi 

describes the Lagrange top on the algebra so(4). 
For q = 2 the perturbed Hamiltonian (|2.2(J|) 

H2 = \(jl + Jl + 2Jl) - Vl . 

describes the Kowalevski top on the Lie algebra so(4). As for the usual Kowalevski top [I], 
we can generalize constraints ()3.3j) and describe the Kowalevski gyrostat on so(4) 

The Lax matrices for the Lagrange and Kowalevski tops on so(4) were constructed in 0j. 
In this section we calculate the corresponding r-matrix ()3.14j) and propose generalization on 
the case p 7^ 3. 

4 Summary 

We construct the Lax matrices for an integrable deformation of the known integrable system 
of two interacting p-dimensional and g-dimensional spherical tops. One gets r-matrices 
associated with this model and with some reduced systems such as generalized Lagrange 
and Kowalevski tops. 
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The matrix y (12. 9|) . which determines perturbation of the known Lax matrix, is 

the solution of the equations (j2.10j) . Relation of such equations on the associative algebras 
with a general theory of the classical /^-matrices on the twisted loop algebras remains an 
open question. Some another constructions of such deformations may be found in [2]. 

The author is grateful to M.A. Semenov-Tian-Shansky , V.V. Sokolov and I.V. Komarov 
for useful discussions. The research was partially supported by RFBR grants 02-01-00888 . 
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